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ABSTRACT 

The standard model of Gamma-Ray Bursts afterglows is based on 
synchrotron radiation from a blast wave produced when the relativistic ejecta 
encounters the surrounding medium. We reanalyze the refreshed shock scenario, 
in which slower material catches up with the decelerating ejecta and reenergizes 
it. This energization can be done either continuously or in discrete episodes. 
We show that such scenario has two important implications. First there is 
an additional component coming from the reverse shock that goes into the 
energizing ejecta. This persists for as long as the re-energization itself, which 
could extend for up to days or longer. We find that during this time the overall 
spectral peak is found at the characteristic frequency of the reverse shock. 
Second, if the injection is continuous, the dynamics will be different from that 
in constant energy evolution, and will cause a slower decline of the observed 
fluxes. A simple test of the continuously refreshed scenario is that it predicts a 
spectral maximum in the far IR or mm range after a few days. 

Subject headings: Gamma-rays: Bursts - Hydrodynamics: Shocks - Radiation 
Mechanisms 



1. Introduction 

The standard model for GRB afterglows assumes that relativistic material is 
decelerating due to interaction with the surrounding medium. A shock wave is formed 
heating the surrounding matter to relativistic temperatures. It is assumed that both 
magnetic fields and accelerated electrons aqcuire an energy density which is a significant 
fraction of the equipartition value. In the simplest case, which will be referred to as the 
standard scenario, a single value of the energy and the bulk Lorentz factor is injected either 
as delta or as a top-hat function, of duration short respect to the afterglow. The total 
energy is fixed in time and equals the initial energy of the explosion. 
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Slower moving material is essential to all models which use density gradients as a means 
of acceleration. Actually, if this is indeed the mechanism, most of the system's energy is 
carried by the slower material. This scenario overcomes the need for a clean environment. 
The fact that one now needs a substantially higher energy input can be addressed by a 
very energetic source such as a massive star. A similar situation exists in some cases of 
supernova. When a shock wave propagates through the envelope of the star arrives at the 
edge it accelerates, and higher and higher velocities are being imparted to a smaller fraction 
of the mass. 

In GRB afterglows, if such slower material with significant energy is ejected, it will 
affect the evolution in two major ways: first, the system becomes more energetic as time 
passes (refreshed shock scenario), therefore the temporal decay of the afterglow will be 
slower (Rees & Meszaros 1998). Second, since the reverse shock will last for as long as the 
energy supply continues, it adds an additional long-living (reverse) emission component, 
typically at low frequencies. The emission from such a reverse shock was considered by 
[Kumar fc Piran 1999| for the discrete injection case. With accurate enough observations 
of the afterglow temporal decay or good spectral sampling, especially at radio to mm 
frequencies, both of these features may be detected, and could therefore constrain the 
possibility of additional energy injection. 



2. Dynamics 

We assume here that the source ejects a range of Lorentz factors. The simplest 
description is that there is a certain amount of mass M(> 7) oc 7~ s moving with a Lorentz 
factor greater than 7, all ejected at essentially the same time ( i.e., over a period which 
is much shorter than the afterglow timescales). The energy associated with that mass is 
7MC 2 oc 7 _s+1 . Note that this is valid only for s > 1 , were for s < 1 the energy above 
any Lorentz factor is constant since it is all concentrated near the highest Lorentz factor. 
We normalize this proportionality using the initial Lorentz factor, 70, and the initial energy 
content, which is about the "burst" energy Eq\ 

£(> 7) ~ £o(7/7o)- s+ \ (1) 

down to some value 7 m i„ < j - This tilted top-hat injection leads to a "refreshed" shock 
scenario ( |Rees fc Meszaros 1998| ), in contrast to the standard straight top-hat or delta 



function model with a mono-energetic E a and a single value 7 Q . Although 70 and E are free 
parameters, we have lower limits on 70 > 100 to avoid pair creation during the GRB itself, 
and we have an estimate of the initial energy Eq from the energy seen in the burst. The 
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actual value of 70 can be obtained from the onset time of the afterglow or from the reverse 
shock initial frequency. This scenario is related to, but not identical, to the [Fenimore & 



Ramirez- Ruiz 19§P| model in which a wind with near-steady 7 impacts on a decelerating 



'wall' produced by the outermost shells of material which first make contact with the 
exterior gas. The refreshed shock or tilted top-hat scenario envisages an 'engine' duration 
which is instantaneous compared to the deceleration time, whereas the latter scenario 
assumes that the wind duration is longer than the initial deceleration time. However if the 
'engine' produced a wind whose Lorentz factor decreased with time, the net effect could be 
rather similar to the refreshed shock. 

The low Lorentz factor mass will catch up with the high Lorentz factor mass only 
when the latter has decelerated to a comparable Lorentz factor. At that time the shocked 
material (both reverse and forward shock) has a Lorentz factor 7 satisfying 

£o(7/7or s+1 ~ £ ~ 7 2 #V- (2) 

If we assume that the outer density is p oc R~ 9 , then we have 

7 = 7 o(i?/i?o)- (3 ^ )/(1+s) (3) 

Using the relation t ~ R/rfc we get 

R = R (t/t ) {1+s)/{7+s - 29) (4) 
7 = 7o(tAo)- (3 - 9)/(7+s - 29) , 

where Rq and t are the deceleration radius and deceleration time of the initial material 
with Eq, 70 Note that these scalings differ from those of Rees and Meszaros (1998), who 
obtained s + 1 where we have s in the scalings of expressions after our equation @. For 
s = 1 the energy in the slow material adds up only logarithmically and therefore the above 
expression degenerates to the usual ones of instantaneous injection. However, for s > 1 
where more energy is stored in a slowly moving material, we have a slower decay as more 
and more energy is added to the system as time evolves. Note also that these expression 
are only valid for g < 3. Otherwise, the shocks is accelerating (see Blandford & Mckee for 
3 < g < 4 and Best & Sari for g > 4). However, the most useful values of g are probably 
the constant density case g = and the wind case g = 2. 

Since the additional slow shells catch up with the shocked material once these are 
with comparable Lorentz factors, the reverse shock is always mildly relativistic (Rees & 
Meszaros 1998, Kumar and Piran 1999). The thermal Lorentz factor of the electrons is 
therefore roughly given by the ratio of proton to electron mass Ye ~ e e m p /m e . In the 
forward shock, the thermal Lorentz factor of the electrons is 7/ ~ e e ^m p /m e . 
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The pressure behind the reverse shock is proportional to the density behind the reverse 
shock which is comparable to the density in front of it (a mildly relativistic shock) therefore 

M G+sg-g 

p r OC n r OC — 5— OC t 7+s-2 9 . 

As a check we can see that this is also the pressure at the forward shock which is 
proportional to 7 2 p. The forward and reverse shock have the same bulk Lorentz factor and 
the same pressure, while the forward shock has a temperature which is higher by a factor 
of 7. 



3. Radiation 

The dynamics above determine the bulk Lorentz factor and the thermal Lorentz factor 
of the electrons as function of time. To estimate the resulting synchrotron radiation, one 
also needs an estimate of the magnetic field. For a given magnetic field, the spectrum 
consists of four power law segments separated by three break frequencies v a , v m and v c 
(Sari, Piran and Narayan 1998, Meszaros , Rees & Wijers 1998). Adopting the standard 
assumption that the magnetic energy density is some fraction 6b of equipartition, i.e., 
proportional to the pressure, we have 

_ 1 6+sg-g 
B OC t 2 7+s-2 9 . 

Since the pressure in the reverse and forward shock is identical, the magnetic field will also 
be the same, if the equipartition parameter e# is the same. The difference between the 
forward and reverse shock lies then in the number of electrons (which is larger by a factor 
of 7 at the reverse shock) and their thermal Lorentz factor (which is smaller by a factor of 
7 at the reverse shock). This result in the following general properties, valid at any given 
moment: 

1) The peak flux of the reverse shock, at any time, is larger by a factor of 7 than that 
of the forward shock: F^ max = jF^ max 

2) The typical frequency of the minimal electron in the reverse shock is smaller by a 
factor of 7 2 : v r m = v f Jl 2 - 

3) The cooling frequency of the reverse and forward shock are equal: v r c = v{ = u c . 

4) At sufficiently early time (typically the first few weeks or months) v r ^ < v r ^ and 
< v c- The sefl absorption frequency of the reverse shock is larger than that of the 

forward shock. It is larger by a factor of 7 s / 5 initially, when both are in fast cooling and by 
a factor of 7 s / 5 if both are in slow cooling. 
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Points one and two above agree with those of Kumar and Piran, calculated for the 
discrete case. We have here generalized the result to include the effect of the cooling 
frequency and self absorption frequency. The combined reverse+forward shock emission can 
therefore be one of three types, evolving in time in the following order: 

A) Both reverse and forward shock are cooling fast: u[ < v r a < v c < v r m < v* m . 

B) Reverse is slow cooling, forward is fast cooling: u[ < v r a < v m < v c < u m . 

C) Both reverse and forward shock are in slow cooling, u[ < v r a < v T m < v m < v c . 

These three spectra are presented in figure |l|. As evident from the figure, the forward 
shock dominates the emission at very low and very high frequencies while the reverse shock 
contributes to a spectral "bump" at intermediate frequencies. The peak flux is that of the 
reverse rather than the forward shock. The combined spectrum is somehwat flatter than the 
usual one (which uses the forward shock only). A good sampling of the spectrum, especially 
at low frequencies, can therefore show the existence or non-existence of such a feature. The 
forward shock alway dominates above v > max(z/^, u c ) by a small factor of / y p ~ 2 . Since the 
value of p is close to 2, the forward shock does not radiate much more than the reverse at 
high frequencies. 



- FIGURE 1 - 



In the case of fast cooling we have ignored the effect of the ordered structure of the 
electron's energy behind the shock ( pranot, Piran fc Sari 2000| ), both for the reverse and 



forward shock. This effect will increase the emission at frequencies below the self absoption 
frequency, v < u a , but will not change the qualitative conclusions of this paper. 

The spectrum displayed in figure [l] is valid at any moment if the energy and momentum 
injection is continuous, but also at the moment of impact in the case of the standard top-hat 
injection. However, in the latter case the reverse shock component will rapidly disappear 
as discussed by Sari & Piran 1999 and Meszaros & Rees 1999. If the injection is discrete, 
the dynamics of the forward shock right after the collision will not be affected, and it will 
evolve as in the standard non-refreshed scenario. . 

For the continuous case, the time dependence t~ q of the various quantities is given in 
Table 1, for arbitrary parameters s and g, assuming a spectral shape cx . Above the 
peak Umax = rnin[i/ m , u c ] where the flux has the value F U)Tnax the dependence F v oc t~ a v~^ is 
calculated separately for the slow and fast cooling regimes. 



- TABLE 1 - 



- 6- 



To give more specific numerical examples we specialize to the constant density case, 
where g — 0. We then have 

vl = 2.0 x 10-Hz (1 + ^\t,_^El\^\Um (5) 

z 

s m = 9.i x io-hz (i + ^-^V^WC^r^ (6) 

to 

u c = 2.7 x 10-Hz (1 + z)-V^!iE£^H^)-%3 (7) 

= 2.6mJy (1 + z)e l ^ 2 E, 2 n l J 2 D L %{^ ) ^ ( 8 ) 

£o 

^ = 12mJy (1 + ^-/SeV^^/^f ^ * } ^ (Q) 

For slow cooling, z/ c > v m > u a , the spectral peak is at u max = u m , and synchrotron 
self-absorption occurs at 

v{ = 3.6GHz (1 + zyh^ef^EU^Tih 1 ^ (">) 

to 

= 43GHz (1 + ^-^leXC^Of) 8 ^ (») 

r 

while for fast cooling, z/ m > z/ c > u sa , the spectral peak is at v max = u c and we have 

4 = o.3ghz (i + zy^tM 1 °-o 1/10 ti /2 (f ) SfTl ( 12 ) 

i£ = 0.8GHz (1 + zr^eZXrnf^l T°( f )^ (13) 

Using the normalization of the peak flux and the break points, the flux can be 
calculated at any frequency. Similar to the standard case, it is possible to test the model 
by comparing the temporal decay and spectral slopes. 

In the standard case (mono-energetic instantaneous injection), the flux above this 
frequency is falling with time, while the flux below this frequency is rising with time. The 
additional energy in the varying injection case tends to flatten the decay rate, and for high 
enough values of s can even make it grow. (From equations (1) and (4) one sees that the 
t/to factors in equations (5)-(13) are equivalent to a power of the ratio of the injected to 
initial energy E/E ). Stated differently, one would need a steeper spectral index to give 
rise to the same observed temporal decay in the refreshed scenario. Table 2 summarizes 
for g = the values of the spectral index f3 that can be inffered from a measured temporal 
decay index a in the instantenous and refreshed scenarios, for reverse and forward shocks. 
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- TABLE 2 - 

It can be seen that the spectal indices that need to exmplain a t~ l decay that is 
observed in many bursts are considerably steeper. Some confusion can occur between a 
forward moderately refreshed (s = 2) shock in the slow cooling regime and a fast cooling 
forward shock in the non-refreshed scenario, as these two scenarios predict similar relation 
between a and (3 for nominal values. However, most other regimes are considerably different 
from the standard instanteneous forward shock prediction even if one only has moderately 
accurate spectral data information. 



4. Specific Bursts 

GRB 970508. - In the case of GRB 970508, the observations show a steep increase in the 
optical and X-ray fluxes between one and two days. This can be interpreted in terms of a 
varying injection event, e.g. Panaitescu, Meszaros & Rees 1998, who consider the radiation 
of the forward shock assuming a large value of s, a value of ■y m in — H < 7 < 7o an d a final 
energy Ef = 3E a . After this time the energy is constant, and the afterglow can be fitted 
with a standard mono-energetic afterglow, e.g. Wijers & Galama 1999. How could we test 
the hypothesis that the "jump" between half a day to two days is indeed due to varying 
injection? It turns out that a delayed energy injection (if is the same for the reverse and 
forward shocks) has a very strong prediction. We can estimate the forward shock break 
frequencies at t — 2 days by extrapolating them back from those at day 12, where the 
spectrum is well studied (e.g. |Wijers fc Galama 1999| , |Granot, Piran fc Sari 1999b]) . This 



results in = 1.3 x 10 12 Hz, F u rnax ~ 1.7mJy, u[ = 3 GHz and 7 = 4. Therefore, according 
to the lower frame of figure [I] the reverse shock should have v r m = 80GHz and v r a = 30GHz. 
The reverse shock signature should be the largest between these two frequencies, where 
the flux should exceed by an order of magnitude the simple extrapolation of the forward 
shock model back to day two. At a more observationally accessible frequency of 20GHz 
and 100GHz the flux increase due to the reverse shock should be a factor ~ 5 and ~ 8 
respectively, i.e. a flux of 2.2mJy and 6mJy, respectively. Unfortunately, the observations 
around 2 days at these frequencies where short and therefore of low sensitivity. The 3cx 
upper limit of 6mJy obtained by BIMA is consistent with energy injection. This prediction 
of an additional low frequency component in GRB 970508 is similar to that of [Kumar &| 
Piran 199~9|. However, taking the self absorption frequency into account we have shown 



here that this would not apply to the usually observed radio frequencies 1.4GHz-8.4GHz 
but to the range of 30GHz-80GHz. It is therefore of great value to obtain low frequency 
observations or strong upper limits if such jumps in the optical and X-ray flux will be seen 
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again in the future. If the reverse shock signature is not seen in other comparable bursts 
where a light curve peak is detected at a certain time (1.5 days in this case), a varying 
injection episode can be ruled out as an explanation for this peak (or else, it would imply 

£B,rev <C €B,for)- 

GRB 990123. - In GRB 990123, a bright 9th magnitude prompt flash was seen about 60s 
after the trigger time ( |Akerlof, et q/.19"99|) , attributable to a reverse shock (Sari fc Piran 
|1999a| , [Meszaros fc Rees 1999| with an initial temporal decay F oc t~ 2 , steeper than the 
subsequent decay F oc t -11 which is attributable to the forward shock. Here, an impulsive 
mono-energetic (standard) injection fits better than a tilted top-hat, since e.g a varying 
s = 2,g = injection would predict a decay a = 1, while an impulsive event gives a = 2, 
in good agreement with observations for a spectral slope /3 = 1 in the fast cooling regime 
(the values are 2/3 and 9/8 in the adiabatic regime). Though no spectral information 
was available for this burst in the first hours, a slope (3 <1 can be assumed based on our 
accumulated knowledge from previous afterglows. 



5. Discussion 

The dynamics and emission of the forward and reverse shocks is controlled by several 
factors, including the continuity and nature of the energy and mass input, the possible 
existence of external density gradients, and the strength of the magnetic fields in these 
regions. A contineous injection of energy with a lower Lorentz factor has as its main 
consequence that it tends to flatten the decay slopes of the afterglow, after it has gone 
through the maximum. An external density gradient (e.g. as in a wind with p ext oc r~ s , 
where g ~ 2 might be typical) has the property of steepening the decay. This applies in 
general both to the forward and the reverse shocks. 

We have suggested ways in which one can attempt to discriminate between the 
standard straight top-hat injection of energy and momentum with a single 7 D and E , 
which then remains constant throughout the afterglow phase, and a refreshed scenario, 
where the injection is also brief (e.g. comparable to the gamma-ray burst duration and 
therefore instantaneous compared to the afterglow timescale) but in which there is a varying 
distribution of 7 and of energy during that injection, so that matter ejected with low Lorentz 
factor cateches up with the bulk of the flow on long timescales. The afterglow energy then 
increases with time. Under the simple assumption that both the forward and the reverse 
magnetic fields are equal (e_B )7 . = £b,/) a remarkable prediction is that in all regimes (both 
shocks are fast cooling, reverse is slow cooling and forward shock is fast cooling or both are 
slow cooling) the reverse shock spectrum joins seamlessly, or with only a very modest step 
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oc 7 P_2 , onto the forward shock spectrum, extending it to lower frequencies. (This could be 
modified if, for instance, e^ r <C cb/, which would give a spectrum with a more pronounce 
through separating the reverse and forward components). Specifically, in the case of GRB 
970508, if a descrete episode of injection that produced refreshed shocks at t ~ 1.5 days 
is the explanation of the step in the X-ray and optical flux at this time, then one would 
expect (for equal forward and reverse e#) an even more dramatic rise of the 20 GHz and 
100GHz flux at 1.5-2 days. This should be about a factor of 5 and 8, respectively, larger 
than expected from the forward shock values extrapolated back to 2 days. In the case of 
GRB 990123 a reverse shock appears to have been responsible for the prompt optical flash, 
and the decay indices (for reasonable spectral slopes) are compatible with monoenergetic 
impulsive (or straight top-hat) injection. 

For future GRB afterglow observations, the main prediction from having comparable 
values of es in the forward and reverse shocks of baryon loaded fireballs is that the peak 
flux is found at the peak frequency of the reverse, rather than of the forward shock, i.e. at 
lower frequencies than typically considered. The IR, mm and radio fluxes would therefore 
be expected to be significantly larger than for simple (forward shock) standard afterglow 
models (e.g. Fig 1). This holds whether the injection in contineous or descrete. The two 
contributions continue to evolve as a pair of smoothly joined components, the ratio of the 
two peak frequencies v^/^L ^ 7 -2 an d peak fluxes F u ^ mr /F^ m f oc 7 gradually approaching 
each other until they coincide at the transition to the non-relativistic case 7 ~ 1. 
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Table 1: Temporal exponents of the peak frequency u m , the maximum flux F„ m , the cooling 
frequency v c and the flux in a given bandwidth F v , calculated both in the adiabatic regime 
v m < v < v c (F v oc F Vm {y m jv)^ oc t~ Q V~ /3 , where f3 = (p — l)/2), and in the cooling regime 
u c < v m < v [F v oc [v c lv m ) x l 2 [v m lvY oc t~ a v~P where (3 = p/2). 
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Table 2: The spectral indices (3 that would be deduced from an observed temporal decay 
index a. The third column gives the well known "standard model" results for instantenous 
injection for which the reverse shock does not last. The fourth and fifth columns are for 
contineous injection case with s = 2 and s = 3 respectively. For extended emission the 
reverse shock lives as long as the slower material keeps arriving, e.g. days. The value in 
brackets demonstrate the numerical value for a nominal temporal index of a — 1. 
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Fig. 1. — Spectrum of the reverse and forward shocks with synchrotron peaks at v r m and u^, 
synchrotron self-absorption frequencies v r a and v[ and cooling frequency u c (same for both, 
assuming = e r B ), for an electron injection spectrum oc 7 _p . Top: both shocks fast cooling. 
Middle: reverse slow, forward fast cooling. Bottom: both slow cooling. 



